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Abstract. Recently, a metal-insulator transition (MIT) was found in the anisotropic Ander- 
son model of localization by transfer-matrix methods (TMM) jl], |j| [J . This MIT has been 
also investigated by multifractal analysis (MFA) [0] and the same critical disorders W c have 
been obtained within the accuracy of the data. We now employ energy level statistics (ELS) 
to further characterize the MIT. We find a crossover of the nearest-neighbor level spacing 
distribution P(s) from GOE statistics at small disorder indicating metallic behavior to the 
Poisson distribution at large disorder characteristic for localized states. An analysis of the 
system size dependence of the spectral rigidity As(L) confirms the values of W c found in 
Ref. |,|. 
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1 Introduction 

Recent TMM studies ||, |[ [| of the anisotropic Anderson model show that an MIT 
exists even for strong hopping anisotropy 7. The values of the critical disorder in 
the band center were found to follow a power law W c oc (1 — 7)^ independent of the 
orientation of the quasi-lD bar. (3 was argued to be independent of the strength of the 
anisotropy. This is supported by multifractal analysis of the eigenfunctions Q , where 
the system size dependence of the singularity spectra is used to determine W c . In 
order to check this further, we employ energy level statistics as another independent 
method. ELS is based on random matrix theory Q and was successfully used to 
investigate the MIT in the isotropic case [|| 0, |j. Directly at the MIT, the ELS 
was argued to be given by a universal distribution P c (s), which should correspond 
to a "critical ensemble" || distinct from the case of a Gaussian orthogonal ensemble 
(GOE) and also from the Poissonian case. However, recent results |j) show that P c (s) 
depends on the boundary conditions and the shape of the samples considered. 
In this work, we show that the MIT in the anisotropic systems can be characterized 
conveniently by ELS and wc find that the critical disorders W c are in good agreement 
with the results of Refs. ||, Q. We further find that P(s) at the MIT in the anisotropic 
systems depends on the anisotropy 7. Thus an estimation of W c by use of the "critical 
ensemble" is incorrect. 
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2 The anisotropic Anderson model of localization 

The Anderson Hamiltonian is given as fl0|| 



ff = X>|i)(*1 + X>iK><il ■ (i) 

We use a simple cubic lattice of size N 3 with orthonormal states \i) at site i = 
(x,y,z). The potential site energies ej are random numbers, uniformly distributed 
in the interval [— W/2, +W/2}. The transfer integrals are restricted to nearest 
neighbors. They depend only on the spatial direction, thus tij — t x , t y or t z . We 
study the two cases of: (i) weakly coupled planes with t x = t v = 1, t z = 1 — 7 and 
(ii) weakly coupled chains with t x = t y = 1 — 7, t z = 1. The anisotropy parameter 7 
ranges from 7 = 0, the isotropic case, to 7 = 1 where the planes/chains are completely 
uncoupled. 

We use the Lanczos algorithm [jfl] to compute the spectrum of H. It is well suited for 
the diagonalization of our sparse matrices |l2| and allows us to compute all eigenvalues 
of H for system sizes of N — 48 on a parallel machine within 60 hours. We use 
50% of the eigenvalues around the band center E = and average over up to 400 
realizations of the random potential, such that at least 2 x 10 5 eigenvalues are used 
for each set of parameters {W,^,N}. Due to the large computational effort, we 
restrict the systematic investigations to sizes up to N = 30. For comparison with 
predictions of random matrix theory, we unfold the spectra by fitting cubic splines 
to the integrated density of states. This sets the mean-level spacing to one. We 
then characterize the local spectral fluctuations by means of the nearest neighbor level 
spacing distribution P(s) and the A3 statistics. The latter measures the rigidity of 
the spectra ||]. 

3 Results 

Extended states in a metal contribute to charge transport even at T = 0. The overlap 
of the extended states results in level repulsion and their spectral properties are 
characterized by the GOE. On the other hand, localized states cannot contribute to 
charge transfer at T = 0, resulting in insulating behavior. The energy levels are 
uncorrelated, consequently the probability that energy levels are close together is 
very high and the ELS is given by the Poissonian statistics. Thus a change from the 
GOE behavior to Poissonian may indicate the existence of an MIT. 
As expected from the isotropic case, we find P(s) to be close to the GOE statistics at 
small disorder and close to Poisson statistics at large disorder as shown in Fig. |l| for 
an already quite strong anisotropy. For stronger anisotropy 7, the transition occurs 
at smaller values of the disorder parameter W. 

In an infinite system, there is a sharp transition from extended to localized behavior 
at the MIT. However, in any finite system, the characteristic lengths scales of the 
states close to the transition will exceed the system size. Thus for a given TV, one finds 
characteristic deviations which result in a continuous change from GOE to Poissonian 
statistics as W is varied across the MIT. Only directly at W c the statistical properties 
are independent of N, because of the scale invariance of the multifractal wave functions 
at the MIT (b| ||. In order to identify the extended, critical and localized regimes 
and to determine the critical disorder W c properly, we therefore examine the system 
size dependence of the ELS. 
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Fig. 1: P(s) for weakly coupled chains with N — 24 and 7 = 0.9. 

As an example we show in Fig. ||the A3 statistics for weakly coupled planes at 7 = 0.9 
for 4 system sizes N ranging from 13 to 30. For W = 6, we find that upon increasing 
N that there is a clear trend towards the GOE prediction. On the other hand, the 
data for W = 12 tend towards the Poissonian result. At W = 9, the A3 statistics is 
independent of N within the accuracy of our calculation. Thus the critical disorder 
for the present example is W c « 9. In addition to this finite-size dependence, we have 




Fig. 2: A 3 (L) for weakly coupled planes with 7 = 0.9 for N = 13,17,21,30. Larger N 
corresponds to larger symbol size. 

also used finite size scaling analysis for quantities such as a — / Q 3 ° A 3 (L)dL. We have 
constructed scaling functions, which further support the values of W c obtained above 
and confirm the one-parameter scaling hypothesis and thus the existence of the MIT 
in these anisotropic systems. Details will be published elsewhere Jl4| . 
The critical disorders obtained by this analysis agree reasonably well with the results 
from TMM [| and MFA j§. W c decreases with increasing anisotropy with a power 
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law W c — 16.3(1 — 7)^, where (3 = 0.25 for weakly coupled planes and (3 ss 0.6 for 
weakly coupled chains, respectively. 

For a given anisotropy and system size, we can also identify a disorder W' at which 
P(s) agrees with the so called "critical statistics" , characteristic for the MIT of the 
isotropic system (cp. Fig. [I]). However, this disorder is much smaller than W c for 
strong anisotropics. And, of course, the value of W changes when we change the 
system size. Thus, P c (s) is not characteristic for the MIT in anisotropic systems. We 
find that upon increasing the anisotropy that the statistical properties at the MIT 
drift slowly from the "critical statistics" P c (s) of the isotropic case || towards Poisson 
statistics. In that sense, the states at W c seem to be less extended in the anisotropic 
system compared to the isotropic case. This coincides with the results of the MFA 
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Fig. 3: W c ("f) for weakly coupled planes (above) and chains (below) as obtained with 
various methods. The thick solid line is the TMM result of Ref. M . The open squares are 
the disorder values, where the "critical statistics" is found for system size N = 21. 
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4 Conclusions 



We find that a metal to insulator transition exists in the anisotropic Anderson model. 
The critical disorders obtained from the ELS coincides reasonably well with the results 
from TMM f| and MFA j§. The system-size independent P(s) and A 3 (L) at the MIT 
depend on the specific values chosen for the microscopic hopping elements t x , t y , and 
t z . They are different for each 7 and the two anisotropy realizations, namely, weakly- 
coupled planes and chains. Furthermore, the ELS at the MIT is also different from 
the ELS of the isotropic case. Thus we find that P(s) at the MIT is not universal, i.e., 
not independent of the microscopic parameters of the model. 

We gratefully acknowledge support by the DFG within the Sonderforschungsbereich 393. 
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